On the ratchet effect 



Ricardo Chacon* 
Departamento de Electrdnica e Ingenierta Electromecdnica, 
Escuela de Ingenierias Industrials, Universidad de Extremadura, 
Apartado Postal 382, E-06071 Badajoz, Spain 

Niurka R. Quintero 
Departamento de Fisica Aplicada I, E. U. P., 
Universidad de Sevilla, Virgen de Africa 7, 4-1011 Sevilla, Spain, 
and Instituto Carlos I de Fisica Teorica y Computacional, 
Universidad de Granada, 18071 Granada, Spain 
(Dated: February 2, 2008) 

PACS numbers: 

Understanding the ratchet effect induced by symmetry breaking of temporal 
forces is a fundamental issue that has remained unresolved for decades. While 
the dependence of the directed transport on each of the ratchet-controlling pa- 
rameters has been individually investigated experimentally, theoretically, and 
numerically, there is still no general criterion to apply to the whole set of these 
parameters to optimally control directed transport in general systems. We re- 
port that, to optimally enhance directed transport, there exists a universal force 
waveform which can be understood as the result of two competing fundamental 
mechanisms: the increase of the degree of breaking of the aforementioned tem- 
poral symmetries and the decrease of the force impulse. We demonstrate that 
this universal waveform explains all the previous experimental, theoretical, and 
numerical results for a great diversity of systems, including motor enzymes. 

Consider a general system (classical or quantum, dissipative or non-dissipative, one- 
or multi-dimensional, noisy or noiseless) where a ratchet effect {1-3) is induced by solely 
violating temporal symmetries of a T-periodic zero-mean ac force / (t) which drives the 
system {4.- 18). This effect is important because of its clear applicability to such diverse 
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fields as biology, and micro- and nano-technology (2,19-23). In particular, it applies to 
electronic transport through molecules (12), smoothing surfaces (24), controlling vortex 
density in superconductors (22), separating particles (25,26), controlling directed current 
in semiconductors (6,9), and rectifying the phase across a SQUID (27). A popular choice 
would be the simple case of a biharmonic force, 



where har\^ represents indistinctly sin or cos, and puj\ = qoo2, P,Q coprime 
integers. In this case, the aforementioned symmetries are the shift symmetry 
(S s : / (t) = —f (t + T/2) , T = 2nq/uJi = 2-Kpjuj 2 ) and the time-reversal symmetries (St r ,± '■ 
f(—t) = ±f(t)). Now a general unsolved problem is to find the regions of the parameter 
space (ej,</?j),ei + e 2 = const., where the ratchet effect is optimal in the sense that the 
average of relevant observables is maximal, the remaining parameters being held constant. 
We show in the following that such regions are those where the effective degree of symmetry 
breaking is maximal. Without loss of generality, we shall illustrate the degree of symmetry 
breaking (DSB) mechanism by using the following working model for the driving force: 



where cn (•; m) and sn (•; m) are Jacobian elliptic functions of parameter m, fl = 2K(m)/T, 
= K(m)9j r n, K(m) is the complete elliptic integral of the first kind, T is the period of the 
force, and 9 is the (normalized) initial phase (9 e [0, 2ir]). Fixing e, T, and 9, the force wave- 
form changes as the shape parameter m varies from to 1 (see Fig. 1). Physically, the moti- 
vation of the choice represented by Eq. 2 is that f e in P (t; T, m = 0, 9) = t sin (2nt/T + 9) /2, 
and that f e u ip (t;T,m = 1,9) vanishes except on a set of instants that has Lebesgue mea- 
sure zero, i.e., in these two limits directed transport is not possible, while it is expected for 
< m < 1. Thus, one may expect in general the average of any relevant observable 3? to 
exhibit an extremum at a certain critical value m = m c as the shape parameter m is varied, 
the remaining parameters being held constant. In this work, we demonstrate that such a 
value m c is universal, i.e., there exists a universal force waveform which optimally enhances 
the ratchet effect in any system. Clearly, there are two competing fundamental mechanisms 
which allow one to understand the appearance of such extremum: the increase of the degree 



f(t) = exharx (uj x t + y^) + e 2 har 2 (oJ 2 t + y? 2 ) , 



(1) 



feiu P (t) = ef(t;T,m,9) 

= e sn (fit + 0; m) cn (fit + 0; m) , 



(2) 
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of breaking of the shift symmetry as m is increased, which increases the absolute value of 
the average, and the effective narrowing of the force pulse as m is increased, which decreases 
the absolute value of the average. The former mechanism arises from the fact that a broken 
symmetry is a structurally stable situation (Curie's principle) and can be quantitatively 
characterized by noting that 

where dn(-;m) is the Jacobian elliptic function. Equation 3 indicates that the degree of 
deviation from the shift symmetry condition (D (t;T,m,9) = 1) increases on average as 
m — > 1, irrespective of the values of the period and initial phase (see Fig. 2). Thus, 
while increasing the shape parameter from m (0 < m < m c ) improves the directed transport 
yielding a higher average, it simultaneously narrows the force pulse lowering the driving 
effectiveness of the force. Indeed, the latter effect becomes dominant for sufficiently narrow 
pulses (m > m c ). Next, one exploits the universality of the waveform corresponding to m c 
to deduce the optimal values of the parameters of the biharmonic force of Eq. 1. To this 
end, notice that we chose the function of Eq. 2 to satisfy the requirement that m c be 
sufficiently far from 1 so that the elliptic force is effectively approximated by its first two 
harmonics. One thus obtains a relationship between the amplitudes of the two harmonics 
in parametric form: ei )2 = £1,2 (j^). This relationship does not depend on the initial phase 
9, and hence neither does it depend on the breaking of time-reversal symmetries of the 
biharmonic approximation corresponding to the elliptic force. For a general biharmonic 
force (Eq. 1), this means according to the DSB mechanism that the optimal ratchet-inducing 
values of the initial phases ip 1 , ip 2 should be those giving a maximal breaking of one of the 
two time-reversal symmetries of the force, while the relationship e 2 = ^2(^1', P, q) should 
control solely the degree of breaking of the shift symmetry. Note that this symmetry is 
not broken when p, q are both odd integers. Consequently, if the DSB mechanism is right, 
the relationship e 2 = ^2(^1', P,q) {p + Q = 2n+ l,n = 1,2,...) controlling the degree of 
breaking of the shift symmetry should be independent of the particular system in which 
directed transport is induced. This implies that any averaged observable < 3? > should 
be proportional to a certain function g (61,62) = g (ei, e^p, q) which is ~ p\ {e\)pi(e-i) in 
leading order with pi(ei) ~ £1,^2(^2) ~ ^ r i s positive integers. Since the aforementioned 
extremum m c is independent of the driving amplitude, one defines e 1 = e (1 — a) , e 2 = ea 
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(a E [0, 1]), so that g (ei, e 2 ) ~ (1 — a) r a s taking e = 1 without loss of generality. Since the 
extremum m c is also independent of the driving period, one has the symmetry relationship 
g (ei, e 2 ;p, g) = 9 (^2, £i; The problem thus reduces to finding the relationship between 
(r, s) and (p,q). From Maclaurin's series, one sees that the simplest function satisfying all 
these requirements in leading order is (1 — a) p a q , and hence g (ei,e 2 ;p, g) ~ efe 2 . Indeed, 
previous theoretical analyses of every kind on a great diversity of systems (5-11,14,15,17,18) 
have found that the averaged observable is always proportional to such a factor in leading 
order. In particular, one thus obtains 



sech 



rKjl-m) 
K{m) 
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K(m) 



v ' m 3 K 6 (m) v ' 

for the biharmonic approximation corresponding to the elliptic force (i.e., p = 2,q = 1). 
Therefore, the shape function S(m) is a universal function which controls the breaking of 
the shift symmetry in leading order for the resonance (p,q) = (2,1). It presents a single 
maximum at m — 0.960057 ~ m c for which e 2 = e 2 (ei) ~ 0.396504ei (note that e 2 = ei/2 
for m = 0.983417; see Fig. 3). Since the DSB mechanism is scale-independent, the critical 
value m c could well be defined by a purely geometric condition (Occam's razor) which 
takes into account the two aforementioned competing mechanisms (degree of breaking of 
symmetries and transmitted impulse): A(m = m c )/A(m = 0) = $/2, where A(m) = 
Jq^ 2 feiiip(t',T,m,9 = 0)dt and $ = (\/5 + l) /2 is the golden ratio. This gives m c = 
0.9830783.... The waveform corresponding to this value can be very accurately approximated 
by a sawtooth of unit height (h — 1), unit period (A = Ai + A 2 = 1), critical asymmetry 
parameter a c = Ai/A = 1/4, and critical tangents tan# Cjl = h/Xi = 4, tan6* Ci2 — h/X 2 — 4/3 
with # Cj2 ee 4arctan$. 

We now discuss the application of the DSB mechanism to time-reversal symmetries. For 
the sake of clarity, consider first the case with no dissipation. According to the above 
discussion, the optimal ratchet-inducing values of the initial phases y?i,y? 2 should be those 
yielding a maximal interference- type effect of the two harmonic functions in Eq. 1, i.e., in 
general, \qip 2 — P<fi\ c can solely adopt some of the values {0, n/2, ir, 37r/2} depending upon 
the particular harmonic functions har\^ and the particular values p, q. When dissipation is 
not negligible, two effects are expected: The average < 3? > will decrease as dissipation is 
increased, and the critical values of the initial phases will be shifted from their values at the 
limiting dissipationless case (ip ci — > tp c i + ip ). Note that the DBS mechanism implies that 
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the dissipation phase ip — > as dissipation vanishes. 

We found that the present theory confirms and explains all the previous experimen- 
tal, theoretical, and numerical results for a great diversity of systems subjected to a 
biharmonic force (Eq. 1) (5-11,14,15,17,18). In particular, it explains recent experi- 
mental results of directed diffusion in a symmetric optical lattice (13), where the force 
F = (constant) [(1 — B) cos (out) + B cos (2ut — <fi)] yielded the maximum velocity of the 
centre-of-mass of the atomic cloud at B ~ 0.33, = ir/2, and the complete data se- 
ries fitted the functional form (1 — B) 2 B (i.e., efe 2 ), in confirmation of the predictions 
above. The analysis of the problem of electron tunneling through a one-dimensional po- 
tential barrier in the presence of an externally applied ac force fi sin (out) + f 2 sin (2uot + tp) 
(11) showed that the directed current is proportional to 2Aj\f 2 cos tp, in confirmation of the 
present theory. Also, the DSB mechanism works powerfully in the case of electron trans- 
port through a molecular wire weakly coupled to two leads and subjected to a biharmonic 
(A 1 sin (Qt) + A 2 sin (2£lt + 0)) laser field (12). Quantum calculations showed (12) that 
only when the symmetry breaking is maximal (i.e., A 2 = Ai/2,(p = n/2) does the average 
current reach its maximum value, and that the average current is proportional to the cou- 
pling strength. For weak coupling, this maximum is about five orders of magnitude higher 
than the corresponding value at = 0. We also found that the universality of the present 
theory is confirmed by previous studies on systems subjected to sawtooth-shaped excitations 
(26,28). Indeed, a study of horizontal transport of granular layers on a vertically vibrated 
sawtooth-shaped base (26) showed that the horizontal transport (v x ) versus asymmetry pa- 
rameter a has an extremum at a = 0.25 only when the sawtooth waveform approximated 
the above universal waveform (i.e., A = 6 mm, h = 5 mm, cf. (26); note that the universal 
waveform is recovered for h = 6 mm). Although we have here limited ourselves to tempo- 
ral symmetries, the DSB mechanism applies to spatial symmetries as well. In this regard, 
one should expect that the ratchet potential underlying biological motor proteins might be 
optimized according to the DSB mechanism. A clear experimental confirmation of this pre- 
diction appears in the context of the actomyosin dynamics in the presence of external loads 
(28), where the mean first-exit time versus external applied force was studied for a sawtooth 
potential profile with different values of the asymmetry parameter. The experimental data 
were optimally fitted only when the sawtooth waveform approximated the above universal 
waveform (i.e., tan#i = h/Xi ~ 4.55, cf. (28); recall that tan# C) i = 4 for the universal 
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waveform) . 

We here studied the example of kink-asisted directed energy transport in a driven, 
damped, sine-Gordon (sG) equation (see Fig. 4), where directed energy transport is pre- 
dicted for a non-zero topological charge, implying the existence of sG solitons (kinks) in the 
system. As Fig. 5 shows, the theoretical predictions are confirmed by numerical simula- 
tions even in the presence of noise. We should stress that the consequences of the present 
theory extend beyond the problem of directed transport. It applies, for example, to the 
phenomenon of synchronization of arrays of coupled limit-cycle oscillators (29), where the 
maximal symmetry breaking of a homogeneous, time-delayed, periodic coupling gave the 
maximum decrease of the synchronization frequency (see Fig. 6). Also, the DBS mech- 
anisms explains the effectiveness of harmonic mixing signals (with frequencies w, 2oo and 
amplitudes r]^,^^, respectively) in controlling (enhancing and suppressing) stochastic res- 
onance phenomena (30,31). In particular, it explains the dependence of the output power 
at the signal frequency uo (P 2 ~ vl, r l2u^ c f- (30)) in a modified Schmitt trigger electronic 
circuit, and the dependence of the time- and noise- averaged mean value of the response 
(To ~ r it r l2ui c f- (31)) of an overdamped two-well Duffing oscillator. We should stress that 
the present theory can be directly applied to optimizing the effectiveness of a two-colour 
laser field for strong harmonic generation [32] as well as to optimally designing synthetic 
molecular motors [33]. 

* Author to whom all correspondence should be addressed. 
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Figure Captions 

Figure 1. Function f(t;T,m,9) (Eq. 2) vs t/T for 9 = and three shape parameter 
values, to = 0, 1 — 1CT 6 (cyan), and 0.96 (magenta), showing an increasing symmetry- 
breaking sequence as the pulse narrows, i.e., asm^ 1. 

Figure 2. Deviation from the shift symmetry condition D (t; T, to, 9 — 0) = 1 (Eq. 3) 
showing an increasing deviation asm^ 1. 

Figure 3. Universal shape function S(m) in leading order for the resonance (p, q) = (2, 1) 
(Eq. 4) exhibiting a maximum at to — 0.960057. 

Figure 4. Average velocity of the kink centre-of-mass versus shape parameter for the sG 
equation U tt - U xx + sin (U) = -(3U t + f e iii P (t) with e = 0.1, (5 = 0.05, T = 20vr, and 9 = 0. 
Previous results (14) from a collective coordinate approach with two degrees of freedom, 
X(t) and l(t) (respectively, position and width of the kink), can be directly applied to obtain 
an ODE governing the dynamics of these two collective coordinates in the presence of the 
elliptic force (Eq. 2): P = -f3P-qf elHp (t), I = I /(2t) + \/{2al) -(51 - (fi|//2)(l + M " 2 P 2 ), 
where the momentum P(t) = M l X /l(t), = y/l2/ (irl ) is the Rice frequency, a = 7r 2 /12, 
and M = 8, q = 2n, and Iq — 1 are, respectively, the dimensionless kink mass, topological 
charge, and unperturbed width. For the biharmonic approximation corresponding to the 
elliptic force (2), one straightforwardly obtains the following estimate for the average velocity 
of the kink: (x (t)^j = e 3 F(j3,T)S (to), where S(m) is the shape function (Eq. 4) and 
F(f3,T) provides the dependence upon the dissipation and the period (14)- The solid, long- 
dashed, and dot-dashed lines represent, respectively, the results from the numerical solution 
of the sG PDE, the results from the numerical solution of the collective coordinates ODE, 
and the above analytical estimate (^X(t)^. The curves corresponding to the two kinds of 
numerical solution present minima at to = 0.983417 (solid line; recall that to c = 0.9830783) 
and to = 0.97904 (long-dashed line), respectively, while the curve corresponding to the 
analytical approximation has its minimum at to = 0.960057. 

Figure 5. Average velocity of the kink centre-of-mass versus shape parameter for the sG 
equation U tt - U xx + sin (U) = -0U t + feiu P (t) + v^V OM) with e = 0.1, (3 = 0.05, T = 
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20-7T, D = 0.001 and = 0, where rj(x,t) is a Gaussian white noise of zero mean and 
correlations {rj (x, t) rj (a/, t')) = 5 (x — x') 5 (t — t'). The lines represent the noiseless cases 
described in Fig. 4, while circles represent the results from the numerical solution of the 
noisy sG PDE which presents a minimum at m — 0.986525 (recall that m c = 0.9830783). 

Figure 6. Symmetry-breaking-induced frequency suppression versus shape parameter for 
the system of N coupled oscillators d^it) / dt = u + 2e ^ f e m p [(f) j (t — r) — <p i (t) ; T, m, 9\ 
with T = 27r, 9 = 0, and where e is the coupling constant, r is the delay, and u is the 
intrinsic frequency of the oscillators. The lowest stable frequency associated with the syn- 
chronization states (0j = Qt + fio) f° r the biharmonic approximation corresponding to the 
elliptic coupling is given by Q m i n ~ c<Jo/[l+27r 2 nre(yf(m)], where n is the number of neighbours 
and aim) = — ; < sech 7rJ ^, 1 ~" 1 ' 1 +4sech ^ZEMlzHi) L The curve exhibits a minimum 

V ' inK^(m) I K(m) K[m) \ 

at m = 0.9845 (recall that m c = 0.9830783) for n = 4, r = 0.1, and e = 3. Note that this 
minimum is the same for all values of r, n, and e, in confirmation of the present theory. 
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